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[l] One of the major types of solar eruption, coronal mass ejections (CMEs) not only impact space 
weather, but also can have significant societal consequences. CMEs cause intense geomagnetic storms 
and drive fast mode shocks that accelerate charged particles, potentially resulting in enhanced radiation 
levels both in ions and electrons. Human and technological assets in space can be endangered as a result. 
CMEs are also the major contributor to generating large amplitude Geomagnetically Induced Currents 
(GICs), which are a source of concern for power grid safety. Due to their space weather significance, 
forecasting the evolution and impacts of CMEs has become a much desired capability for space weather 
operations worldwide. Based on our operational experience at Space Weather Research Center at NASA 
Goddard Space Flight Center (http://swrc.gsfc.nasa.gov), we present here some of the insights gained 
about accurately predicting CME impacts, particularly in relation to space weather operations. These 
include: 1. The need to maximize information to get an accurate handle of three-dimensional (3-D) CME 
kinetic parameters and therefore improve CME forecast; 2. The potential use of CME simulation results 
for qualitative prediction of regions of space where solar energetic particles (SEPs) may be found; 3. The 
need to include all CMEs occurring within a ~24 h period for a better representation of the CME 
interactions; 4. Various other important parameters in forecasting CME evolution in interplanetary space, 
with special emphasis on the CME propagation direction. It is noted that a future direction for our CME 
forecasting is to employ the ensemble modeling approach. 

Citation: Zheng, Y., et al. (2013), Forecasting propagation and evolution of CMEs in an operational setting: What 
has been learned. Space Weather, 11, 557-574, doi:10.1002/swe.20096. 


1. Introduction 

[ 2 ] The three main types of solar transient/structure 
responsible for space weather are: 1. Solar flares, 2. 
Coronal mass ejections (CMEs), and 3. Corotating interac- 
tion regions (CIRs) and the accompanying high-speed 
solar wind streams. 

[ 3 ] Solar flares and CMEs are arguably the two most 
energetic/catastrophic solar phenomena [e.g., Hudson, 2011; 
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Webb and Howard, 2012; Qiu et al., 2004; Qiu and Yurchyshyn, 
2005; Yashiro et al., 2005; Wang and Zhang, 2007; Zhang et al., 
2001 and references therein]. A CME is the spectacular 
ejection of mass and magnetic field outward from the solar 
atmosphere into interplanetary space. Although it is often 
associated with a solar flare, a one-to-one correlation is not 
always observed. Interplanetary CMEs (ICMEs) are the 
interplanetary manifestations of CMEs typically remote 
sensed by coronagraphs. The speed of CMEs ranges from 
~200km/s to -3000 km/s, with the mass often totaling more 
than 10 12 kg and the kinetic energy exceeding 10 15 J. Often 
flare (electromagnetic) radiation energy can be comparable 
to CME kinetic energy. The relationship between a flare 
and its associated CME has been interpreted in various ways. 
Some studies suggest that rather than being distinct physical 
processes, they could be merely different manifestations of 
the same physical process [e.g., Zhang et al, 2001; 2004], 

[ 4 ] In comparison with solar flares and CIRs, CMEs are 
more likely to have significant (at times severe) space 
weather impacts. Earth-directed CMEs (often originating 
from the disk center — within a longitude range of ±30°) 
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are the drivers of major geomagnetic storms [e.g.. Cane 
and Richardson, 2003; Gopalswamy et al., 20005a, 2005b; 
Gopalswamy, 2006], capable of causing rapid and large 
magnetic field disturbances, enhanced currents in space 
and on the ground (the so-called Geomagnetically 
Induced Currents — GICs), and enhanced radiation in 
our near-Earth environment — both from energetic 
electrons and from ions. CMEs pose serious threats to 
spacecraft components and operations, communication 
and navigation systems, astronauts either in extravehicular 
activity or behind a protective shielding [Lanzerotti, 2004], 
as well as ground systems such as power grids and 
pipelines [e.g., Boteler et al., 1998; Pirjola, 2007]. 

[ 5 ] Another aspect of CMEs' importance in space 
weather lies in their roles in accelerating solar energetic 
particles (SEPs). Large SEP events ("large" means high 
particle flux for protons (or other ion species) at energies 
equal or above 10 MeV here) constitute a serious radiation 
hazard to spacecraft systems and humans in space. 

[6] Physical understanding of SEPs has evolved signifi- 
cantly over the years. Originally, solar flares were believed 
to be the source of the energetic particles; such a view 
largely reflected a lack of CME observations. A paradigm 
shift with regard to SEPs occurred around the 1980s and 
1990s [e.g.. Realties, 1999]. Attention was now drawn to the 
importance of CME-driven shocks for SEP acceleration: 
the shock accelerates and injects particles continuously as 
it propagates outward. It was then suggested that SEP 
events be classified into one of two categories: impulsive 
or gradual. In gradual events, particles are accelerated by 
the fast-mode CME shocks in the outer corona and in 
interplanetary space, while impulsive events are acceler- 
ated at the flare site/active region in the solar corona. 

[ 7 ] Later observational evidence of the "mixed" events 
exhibiting distinct signatures of both impulsive and grad- 
ual events [e.g., Cohen et al., 1999, Cane et al., 2003, 2006] 
raised questions about the effectiveness of such a dichot- 
omy in classifying the SEP events [Cliver, 2008, Kallenrode, 
2003]. Kallenrode [2003] has suggested that pure impulsive 
and pure gradual events are only the extreme cases of 
SEP events and there is a continuous transition between 
different kinds of SEP events. 

[8] Though the relative roles of the flare and CME in SEP 
acceleration remain an active research topic, mounting 
evidence indicates that the CME-driven shock plays a 
dominant role in generating large SEP events [e.g., Cliver 
and Ling, 2009; Gopalswamy et al., 2003; Tylka and Lee, 
2006]. Being able to model a CME and its angular extent 
accurately will not only aid in assessing its impacts on a 
specific spacecraft or a region of interest, but also aid 
in estimating the spatial domain where an SEP event 
is expected to be observed, as illustrated by Figure 5 in 
Li and Zank [2005]. 

[ 9 ] This paper is organized as follows. Section 2 provides 
a brief description of the model used for modeling the evo- 
lution of CME(s) in interplanetary space — the coupled 
Wang-Sheeley-Arge (WSA)+ENLIL+Cone model. Section 
3 describes what we have learned about accurately 


modeling CMEs and their impacts based on our space 
weather operations. Summary and discussion are presented 
in section 4. 

2. The WSA+ENLIL+CONE Model 

[ 10 ] A number of models are available for simulating the 
propagation and evolution of CMEs in interplanetary 
space, each with its own strengths and weakness. For 
example, the coupled coronal (a resistive MHD model — 
See Linker and Mikic [1995]) and heliospheric (a global heli- 
ospheric model called ENLIL) model suite described in 
Odstrcil et al. [2002] is capable of capturing many large-scale 
features of the solar corona and heliosphere. The model 
result was used to interpret the global context of CME obser- 
vations by ACE at 1 AU and Ulysses at 5 AU [Riley et al., 2003]. 
Other CME models include the EIAF (Hakamada-Akasofu- 
Fry) v2 [e.g., Wu et al., 1983; Fry et al., 2001] + 3-D MHD 
Model Ensemble code as described by Wu et al. [2011, and 
references therein] and the MHD model in Shen et al. [2011] 
and Zhou et al., [2012], the coupled WSA+ENLIL+Cone model 
[e.g., Odstrcil et al., 2004a, 2004b, 2005], the drag-based 
analytical CME propagation model [e.g., Vrsnak et al., 2010], 
the Solar Corona, Eruptive Event Generator, and Inner 
Heliosphere modules of the Space Weather Modeling 
Framework (SWMF) [Tdf7i et al., 2005] and so on. It is not the 
intent of this paper to provide an exhaustive list of all the 
available CME models. The rationale for our model selection 
is discussed below. 

[ 11 ] March 2010 marked the establishment of the Space 
Weather Research Center (SWRC) at NASA Goddard Space 
Flight Center (GSFC), a part of the Community Coordinated 
Modeling Center (CCMC, http://ccmc.gsfc.nasa.gov), with a 
primary objective of addressing the space weather needs of 
NASA's robotic missions through experimental research 
forecasts, notification, analysis, and education (http://swrc. 
gsfc.nasa.gov). The center represents maximal leverage of 
the latest scientific research results and more than a decade 
of modeling capabilities enabled by CCMC. SWRC aims to 
be a locus of innovation in research-based space weather anal- 
ysis and forecasting prototypes. It not only addresses the 
unique space weather needs of NASA users, but also pioneers 
next generation space weather prediction systems. CCMC / 
SWRC differs from NOAA Space Weather Prediction Center 
(SWPC) in many aspects; one important difference is that 
CCMC/SWRC broadens the conventional space weather re- 
gime to include interplanetary space weather [Guhathakurta, 
2013]; i.e., it monitors/tracks CMEs (together with other space 
weather effects) in all directions, while NOAA/SWPC is most 
concerned with Earth-directed eruptions and space weather 
effects in geospace. The disclaimer notice posted on the 
CCMC/SWRC website defines these different roles and 
points out that NOAA/SWPC provides the official U.S. 
Government space weather forecast. 

[ 12 ] Taking into account previous research, model vali- 
dation results, model performance [e.g., Odstrcil et al., 
2004a, 2004b, 2005; Taktakishvili et al., 2009, 2011], and run- 
ning time, our space weather operations use the coupled 
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WSA+ENLIL+Cone model, which is also running at 
NOAA/SWPC [Pizzo et al., 2011]. As the most sophisticated 
model currently available to space weather forecasters, it is 
capable of providing a 1-2 day lead-time forecasting for 
major CMEs. Of course, as with all current models, it also 
has limitations as evidenced by discrepancies when com- 
pared with observations. The discrepancies can be attrib- 
uted to multiple factors, including the idealistic "cone" 
representation of CMEs in the model and the absence of 
an internal magnetic structure in the modeled CMEs. The 
latter is a major limitation of almost all models; adding 
the CMEs' internal structure remains both physically 
challenging and computationally expensive. It should be 
stressed that the focus of this paper is not on fine tuning 
the model inputs to produce the best agreement with 
observations, but rather on what has been learned by 
performing real-time CME runs. 

[ 13 ] The Wang-Sheeley-Arge (WSA) model [Schatten, 
1971; Wang and Sheeley, 1995] is a combined empirical and 
physics-based representation of the quasi-steady global 
solar wind flow that can be used to predict ambient solar 
wind speed and Interplanetary Magnetic Field (IMF) 
polarity at Earth. The coronal field is determined by a 
potential field source-surface calculation, for which the 
inner boundary is specified with the use of synoptic 
magnetograms [Arge and Pizzo, 2000]. WSA computes the 
solar wind speed at the source surface using an empirical 
relationship. Global Oscillation Network Group (GONG) 
solar magnetograms are used as input for the WSA model 
in our operation. 

[ 14 ] ENLIL is a time-dependent, 3-D ideal MHD model 
of the solar wind in the heliosphere [Odstrcil et al., 2002, 
2004a, 2004b], designed to treat supersonic outflows in the 
limit where resistivity and viscosity are minimal. It is based 
on the polytropic equation of state. The ENLIL model 
domain extends from the solar equator to within ±60°. 
This grid has sufficient latitudinal range to minimize 
the effects of neglecting high-latitude behavior on low 
latitudes [Odstrcil and Pizzo, 1999]. 

[ 15 ] The cone model concept (similar to Zhao et al. [2002] 
and Xie et al. [2004]) is used to model a CME and its pro- 
pagation through interplanetary space, assuming that it 
propagates with nearly constant angular width in a radial 
direction and that the expansion is isotropic. Although rather 
simplistic, it has been reasonably successful in representing 
CMEs and their impacts at points of interest in interplanetary 
space [Taktakishvili et al., 2011; Palkenberg et al., 2010; 2011] in 
combination with the WSA+ENLIL model. 

[ 16 ] In the coupled WSA+ENLIL+Cone approach 
[Odstrcil et al., 2004a, 2004b], output from the WSA is used 
to set up the ENLIL inner boundary condition at 21.5 Rs 
(solar radii) so that relatively realistic solar coronal mag- 
netic field and structures are taken into consideration. 
ENLIL, together with the Cone model, accounts for the 
ambient solar wind, transient CME structures, and their 
interaction and propagation in the simulation domain of 
the heliosphere. The modeling results from the combined 
CME model are not only useful for forecasters, but also 


valuable for interpreting scientific results [e.g., Riley et al., 
2003; Odstrcil et al., 2004a, 2004b; Baker et al., 2013; Mewaldt 
et al., 2013]. 

[ 17 ] The main input parameters for starting the WSA 
+ENLIL+Cone model are the kinematic parameters of 
CMEs, which are obtained primarily from the white-light 
coronagraph images provided by SOHO, STEREO A, and 
STEREO B spacecraft. A robust triangulation method, 
evolved from the one described in Pulkkinen et al. [2010] 
and based on the projection matrix formulation, is used to 
determine the kinetic properties of a CME. The method 
and the corresponding triangulation tool were developed 
for quick and easy use in an operational space weather 
forecasting setting. Coronagraph images from at least two 
spacecraft are needed for the triangulation method. 
During a specified interval of interest, at least two images 
from each of the two spacecraft are needed to derive the 
CME parameters. The leading edge of a CME image is 
used to derive CME speed. The details of the method will 
be reported elsewhere. 

[is] The input parameters for the Cone model include 
the propagation direction of a CME in terms of longitude 
and latitude in Heliocentric Earth Equatorial (HEEQ) 
coordinates, its speed, its angular width, and the time of 
its arrival at 21.5 Rs (ENLIL's inner boundary). 

3. CME Forecasting in an Operational Setting 

[ 19 ] There are a number of challenges that need to be 
dealt with in a real-time setting as opposed to postevent 
CME analysis. Due to the fact that the coronagraph images 
are all from science missions (SOHO, STEREO A, and 
STEREO B), and there are data downlink limitations, data 
gaps in real-time coronagraph images are to be expected. 
In addition, the (near) real-time coronagraph beacon 
images on STEREO A and B are compressed to minimize 
the data downlink size so an inferior quality of images is 
often encountered in an operational setting, although 
many of these problems can be alleviated when dealing 
with science quality data for research purposes. The li- 
mitations and constraints of real-time observations and 
the desire to achieve accurate forecasting of CMEs necessi- 
tate exploiting and synthesizing all information sources to 
better assess CME characteristics. The example below 
illustrates the importance of utilizing a CME's 3-D proper- 
ties in forecasting CME impacts. 

3.1. The Importance of Synthesizing All 
Available Information 

[ 20 ] To improve the forecast of CME propagation/evolu- 
tion, it is necessary to synthesize all data sources available 
to extract as much information as possible. For example, 
the (near) real-time SDO (Solar Dynamic Observatory) 
(http://sdo.gsfc.nasa.gov) images at high spatial and 
temporal resolution, particularly the ALA data in different 
EUV wavelengths, have been invaluable in our operation, 
allowing us to look at the CMEs in different "lights" 
(different EUV wavelengths in addition to the white-light 
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Table 1. Parameters of the CME That Occurred on 4 October 2011 


CME Starting Time 

Speed km/s 

Direction LON/LAT in HEEQ 

Half-Cone Angle (Degrees) 

Time at 21.5 Rs 

2011-10-04 13:10 UT 

1370 

-155/35 

50 

2011-10-04 15:23 UT 


coronagraph images) and providing valuable additional 
information for improving the CME input parameters 
needed for the global Heliospheric model. Similarly, the 
EUV images from STEREO A and B are often used to aid 
our operation (particularly for farside CMEs), along with 
many other types of observations, such as radio emissions. 
Recent research results have shown that CMEs can change 
direction during their propagation [Byrne et al., 2010]. 
During our operations, several CMEs originating from a 
high-latitude active region were found to deflect toward 
the ecliptic plane. A systematic effort is under way to syn- 
thesize all information sources in order to better determine 
CME input parameters (especially when coronagraph 
images are only available from one spacecraft), and thus 
further improve our forecasts. 

[ 21 ] As described by Cane and Lario [2006], energetic par- 
ticle observations in interplanetary space can be used to 
probe the origin, development, and structure of coronal 
mass ejections. The example below shows how we used 
the concurrent SEP observations (in combination with 
others) to assist our CME impact forecast. 

[ 22 ] Table 1 lists the input parameters for the CME that 
occurred on 4 October 2011. This is a backsided, fairly 
high-latitude (35°) CME with a fast speed (-1370 km/s) 
and a broad angular width starting around 4 October 
2011 13:10 UT. Enhancement in energetic particle fluxes 
was detected over a broad longitudinal extent as 
evidenced in Figure 1: both STEREO A and STEREO B 
saw flux increases in the 13-100 MeV protons following 
the onset of the CME. The MESSENGER spacecraft in 
the vicinity of Mercury also detected enhancement in 
SEP fluxes. Because particles tend to follow magnetic 
field lines, they can be used to trace field line topologies. 
Also, shock-accelerated populations provide information 
about the size of a CME shock as it propagates from the 
Sun to the observer. 

[ 23 ] Panels (a) and (b) of Figure 2 provide the simulated 
three-dimensional density map of the solar wind at 7 
October 2011 06:00 UT (interchange with 7 October 2011 
T06:00Z) and 12:00 UT respectively. The solar wind density 
distribution from three different viewing perspectives is 
displayed, with the left panel providing a view of the ecliptic 
plane, the middle panel showing the meridian plane cutting 
through the Sun-Earth line, and the right panel showing the 
Mercator projection of the 1 AU (Astronomical Unit) spher- 
ical surface (the horizontal axis showing the latitude and the 
vertical axis showing the longitude). Our simulation results 
indicate that the CME may barely touch STEREO B (the au- 
tomatic routine indicated a non-encounter with STEREO B). 
In addition, its relatively high-latitude (35°) propagation also 
acts against its potential impact on STEREO B, which lies 
close to the ecliptic plane. However, the observed flux 


enhancement in the 13-100 MeV energy range by STEREO 
B shortly after the onset of the CME suggests strongly that 
there is a good magnetic connectivity of the CME shock 
and STEREO B and the size of the CME shock is broad 
enough to affect STEREO B. By synthesizing all this infor- 
mation (both from model simulations and observations), 
we issued our CME alert (http://l. usa.gov/lflfEjp) and 
provided an estimated arrival time for both MESSENGER 
and STEREO B. In the alert, we stated that "Simulations in- 
dicate that the leading edge of the CME will reach STEREO 
B around 2011-10-07T06:00Z (plus minus 7 hours) - the flank 
of the CME" (Figure 2 indicates that a more accurate time 
would be somewhere in between 06:00 UT and 12:00 UT). 
Our model estimate was close to the actual shock arrival 
around 7 October 2011 11:30 UT, as indicated by the in situ 
plasma and IMF measurements onboard STEREO B (See 
Figure 3, quantities from STEREO B are in blue). 

3.2. Added Benefit of the CME Simulation Results 
for SEPs 

[ 24 ] An active day in terms of solar activities was 7 March 
2011 with occurrences of 5 M-class flares and two signifi- 
cant CMEs. The two significant CMEs are associated with 
the Ml. 9 class flare, which peaked at 14:30 UT (starting at 
13:45 UT, peaking at 14:30 UT and ending at 14:56 UT) from 
the Active Region 11166 (N11E21) and the M3. 7 class flare, 
which peaked at 20:12 UT (starting at 19:43 UT, peaking at 
20:12 UT and ending at 20:58 UT) from the Active 
Region 11164 (N23W50), respectively. Type II and Type IV 
radio emissions were also observed with the Ml. 9 flare at 
14:30 UT, while the Type II radio emission and a Tenflare 
(strong emissions in the 10.7 cm wavelength band) were 
detected in association with the M3.7 flare. 

[ 25 ] Figure 4 shows the real-time SDO EVE (Extreme 
Ultraviolet Variability Experiment) [Woods et al., 2010] 
images of the two flares at times close to their peak. 
Although the quality of these images is not superlative, 
they provide quick visual information on where a flare/ 
CME originates, and are thus useful for space weather 
operations. The location of a flare/CME can help with 
initial assessment of their potential space weather impact 
as past statistical results [e.g., Gopalswamy et al, 2005c] 
indicate that the geoeffective CMEs tend to occur within 
30° of the center meridian, while Earth-impacted SEP 
events tend to favor the vicinity of 60° west (the right side 
on the images) in longitude. The space weather impacts 
of the two flares/CMEs conform to the statistical results 
(more details are provided below), with the first CME 
being more geoeffective and the second CME (and its 
associated flare) being more SEP effective at Earth and is 
unlikely to have strong geomagnetic storm effects. 
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Figure 1. STEREO beacon data of energetic particle flux at different energy channels, with the 
left side ones from STEREO B and the right side from STEREO A. 


[ 26 ] Figure 5 shows the locations of the STEREO space- 
craft (B for STEREO B and A for STEREO A). SOHO is lo- 
cated at the LI Lagrangian upstream location along the 
Sun-Earth line, about l/100 th of the distance from Earth to 
the Sun. Figure 6 contains the coronagraph images for two 
CMEs. The top images are the snapshots of the first CME 
around 7 March 2011 17:00 UT viewed from three different 
spacecraft: STEREO B, SOHO, and STEREO A. The bottom 
row shows the snapshot of the second CME around 7 March 
2011 21:00 UT from three different viewpoints. The time se- 
quence of these coronagraph images from three spacecraft 
was used (together with the ephemeris information of the 


spacecraft) to extract kinetic CME parameters that are 
needed for the WSA+ENLIL+Cone model. 

3.2.1. Nominal Use: Time of Arrival Prediction and 
Impact Assessment 

[ 27 ] The Cone model input parameters for the two CMEs 
are listed in Table 2. In the WSA+ENLIL+CONE model 
suite, CME density can be modified according to each indi- 
vidual CME. But operationally, we use the default setting 
where the CME density is four times of the ambient solar 
wind. The reason is that determining CME density from 
real-time data streams remains difficult in comparison 
with the other parameters mentioned above. This set of 
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Figure 2. (a) The simulated 3-D solar wind density map at 7 October 2011 06:00 UT. (b) The 
simulated 3-D solar wind density map at 7 October 2011 12:00 UT. 


parameters is the same one we used for our space weather 
operations. The simulation results were communicated to 
our users at near real time (normally with a couple of hours 
delay due to the waiting time for enough coronagraph 
images available for analysis and another ~20 min for the 
model execution. The 1-4 days lead-time forecasting 
capability from the WSA+ENLIL+Cone model still makes 


the effort worthwhile). They are not adjusted/fine-tuned 
to match the observational results. 

[ 28 ] The temporal evolution of the solar density, solar 
wind speed from the ENLIL+WSA+Cone modeling results 
(dashed lines) at LI (around Earth) and STEREO A is 
shown in Figure 7, along with the observations (solid lines), 
with Figure 7a showing those at Earth and Figure 7b 
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Figure 3. In situ beacon data of solar wind plasma and 
IMF from STEREO B (blue) and STEREO A (red). Note 
the shock arrival at 11:30 UT detected by STEREO B. 


showing those at the STEREO A location. The vertical 
dashed lines indicate the observed shock arrival time at 
LI or STEREO A, while the dotted vertical lines are the 
corresponding simulated arrival time. 

[ 29 ] It should be reiterated that unlike research studies, 
the simulation results shown here (and throughout this 
paper) are from real-time operations (using the default 
model setting), which do not favor the fine tuning of 
various modeling parameters /settings to match the obser- 
vations. For example, the model does not reproduce the 
first density peak at STEREO A on 8 March, which may 
have potentially slowed down the arrival of the CME at 
STEREO A and explain the earlier simulated arrival time 
in comparison with the observed one. In addition, at LI, 
the shock structure from the simulation is more gradual 
than ACE's observations. 

[ 30 ] From an operational perspective, however, the 
forecast arrival time at Earth (LI) and STEREO A was 
reasonably accurate. ACE observed an interplanetary 
shock at 10 March 2011 06:00 UT (the dashed vertical line 


in Figure 7a), while our predicted arrival is around 10 
March 2011 06:52 UT (the dotted vertical line in Figure 7a). 
Similarly, STEREO A observed an interplanetary shock at 
9 March 2011 06:50 UT as indicated by the IMPACT (In situ 
Measurements of Particles and CME Transients)/PLASTIC 
(PLAsma SuperThermal Ion Composition) beacon data 
(the dashed vertical line in Figure 7b), while our predicted 
arrival is around 9 March 2011 00:45 UT (the dotted vertical 
line in Figure 7b). Both are within the statistical error (± 7 h) 
derived from the geoeffective historical events [Taktakishvili 
et ah, 2009]. Model results, it should be noted, usually lack a 
clear distinction between the shock and the CME material 
itself. The arrival time in our operation is defined when 
the time derivative of the modeled dynamic pressure 
exceeds a certain threshold value. 

[ 31 ] As a part of our CME forecasting, we also estimate 
the maximum value of Kp (a geomagnetic activity index, 
ranging 0-9) if a CME(s) heads toward Earth. Although 
Kp is a rather simple and low temporal cadence index, it 
is used widely in operational space weather. The Kpmax 
is derived based on Newell et ah [2007] formula using the 



Figure 4. SDO EVE images of the (a) M1.9 and (b) M3 .7 
flare at times close to their peak intensity. 


563 



ZHENG ET AL.: FORECASTING CMES 


'■JB i 


0.0 


0.5 


t Wercury 


A. 


iVenus 


t ' 

i ' 


1.0 


I £ort hr ' 


- 1.0 


- 0.5 


0.0 

Y (HEE) 


0.5 


1.0 


Figure 5. Locations of STEREO A and STEREO B in 
HEE (Heliocentric Earth Ecliptic) coordinates in refer- 
ence to the Sun. The SOHO spacecraft is very close to 
Earth in this picture, 1/100 th distance of the total Sun- 
Earth distance. Courtesy of the STEREO orbit tool at 
http://stereo-ssc.nascom.nasa.gov/where/. 


modeled solar and IMF conditions. Since WSA+ENLIL 
+Cone does not give the polarity of the field in a realistic 
manner for a CME(s), we assume variety of IMF clock 
angles. Kpmax results by assuming the maximum of the 
modeled field to be southward. For the run above, the 
estimated Kp is 4 (same as what was observed). 

[ 32 ] While the above example demonstrates the default 
usage of the WSA+ENLIL+Cone simulation results, next 
we will show a rather novel usage of the CME simulation 
results, which provide a qualitative and graphical estimate 
of SEP distributions that have proven to be useful/informa- 
tive for our NASA mission operators. 

3.2.2. CME Simulation Results: Visual Representation 
of SEP Distribution 

[ 33 ] From a space weather perspective, another interest- 
ing aspect of the 7 March event is that an enhancement of 
solar energetic particle (including electrons and ions) fluxes 
was observed across a broad longitude in space, at space- 
craft along the Sun-Earth line by SOHO and ACE (LI), 
and GOES at geosynchronous orbit, and at STEREO A. 
This SEP event is associated with the M3. 7 flare (the flare 
peak time is marked by the left dashed line in Figure 8a) 
and the 1980 km/s CME. Figure 8 shows the temporal profile 
of fluxes at different energy levels from different spacecraft. 


\> > 

♦ 

Jo, 

/ 

201 1/03/07 17 06 

V y 

1 * 

| 4 | JOUO-CMimoc; • | ► 

□ □LiiEiinao m 

fit ii; uahti in ra i » 

'n * 

• ■ 

0 - 

201 1/03/07 21 18 

# 3 

ni cl : j j TH n ^ « » ■ n rrrTrrm n r n 

j (j, *•» 


Figure 6. Coronagraph images for two CMEs. The top images are the snapshot of the first 
CME around 7 March 2011 17:00 UT viewed from three different spacecraft: STEREO B, 
SOHO, and STEREO A. Similarly, the bottom row shows the snapshot of the second CME 
around 7 March 2011 21:00 UT from three different vantage points. 
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Table 2. The Cone Model Input Parameters for the Two CMEs on 7 March 2011 


CME # Associated Flare 

Speed km/s 

Direction LON/LAT in HEEQ 

Half-Cone Angle (Degrees) 

Time at 21.5 Rs 

M1.9 flare 14:30 UT 

710 

-13/15 

35 

2011-03-07 19:51 UT 

M3.7 flare 20:12 UT 

1980 

50/17 

45 

2011-03-07 21:40 UT 


Figure 8a displays those from STEREO A, ACE, and GOES 
in reference to the solar flare intensity in two different X- 
ray bands in units of watts/m 2 (the 1-8 A wavelengths in 
black and the 0.5-4 A in blue at the bottom of the plot). 
Figure 8b shows energetic proton fluxes measured 
by SOHO/COSTEP (Comprehensive SupraThermal and 
Energetic Particle Analyzer). 

[ 34 ] For Figure 8a, the dark red line is the 38-53 keV 
electron differential flux (unit: l/cm 2 /s/sr/MeV) from 
ACE, the green line is the 35-65 keV electron differential 
flux (unit: l/cm 2 /s/sr/MeV) from STEREO A, the purple 
line is the >10MeV proton integral flux from ACE (unit: 
l/cm 2 /s/sr), the red line is the >10MeV proton integral 
flux from GOES 13 (unit: l/cm 2 /s/sr), the black line is the 
2.2-12 MeV proton differential flux from STEREO A (unit: 
l/cm 2 /s/sr/MeV), and the yellow line is the 13-100 MeV 
proton differential flux (unit: l/cm 2 /s/sr/MeV) from 
STEREO A. For Figure 8b, the black line is the 4—9 MeV 
proton differential flux, blue line for the 9-15.8 MeV proton 
flux, light orange for the 15.8-39.8 MeV proton flux, and 
dark orange for the 28.2-50.1 MeV proton flux, all of which 
are in the unit of l/cm 2 /s/sr/MeV. 

[ 35 ] We can see that the SEP event has the characteristics 
of a gradual event in which the CME plays a significant role 
[e.g., Luhmann et ah, 2007; Lion et ah, 2013]. The energetic 
proton fluxes on STEREO A, in particular, show continu- 
ous increases as the CME-induced shock and ejecta move 
out toward STEREO A. When the shock/CME arrived at 
STEREO A, there was a reduction in SEP event intensity 
(indicated by the right dashed vertical line on Figure 8a) 
as if the strong magnetic fields in the CME ejecta form a 
barrier to particles' entry into that structure, often referred 
as the Forbush decrease [e.g.. Cane, 2000]. 

[36] Due to complex physical processes and the different 
makeup (both in composition and energy spectrum) of SEP 
sources (including both flare and CME/ICME shock contri- 
butions) involved in different SEP events, a realistic SEP 
model capable of forecasting/modeling SEP event occur- 
rence in three-dimensional interplanetary space (applica- 
ble for different spatial locations) and at the same time of 
operational use (computationally feasible) remains elusive 
[Luhmann et ah, 2010; Zank et ah, 2007], The simulation 
results from the combined WSA+ENLIL+Cone model can 
be very useful in terms of qualitatively/visually providing 
an estimate of where an SEP event would be expected. 
The activity location plus the simulated CME/shock profile 
(both spatial and temporal) provide information on 
the magnetic connectivity to an observer and therefore 
facilitate a rough estimate of the expected SEP spatial 
extent. This approach assumes that the influence of the 
shock evolution dominates over the diffusive transport in 


determining the SEP profiles as shown by Luhmann et ah 
[2010]. Figure 9 serves as such an example, showing the 
solar wind density distribution from three different view- 
ing perspectives. The “*"s and "x"s indicate where the 
SEPs can be expected based on the simulated CME/shock 
properties and magnetic connectivity (the black and white 
dashed lines are magnetic field lines), which are consistent 
with the SEP observations at STEREO A and along the 
Sun-Earth line very well. 




00:00 00:00 00:00 00:00 00:00 

2011/03/08 03/09 03/10 03/11 

Figure 7. Simulated timeline of solar wind density and 
speed at LI and STEREO A locations (dashed lines) and 
comparisons with observations (solid lines) (Run 1). 
The dashed line in the top panel of Figures 7a and 7b 
represents the density and the dot-dashed line repre- 
sents the speed. The vertical dashed lines indicate the 
observed shock arrival time at LI or STEREO A, while 
the dotted vertical lines are the corresponding simu- 
lated arrival time. 
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Figure 8. (a) SEP measurements from STEREO A, ACE, and GOES 13 in reference to the solar 
flare X-ray fluxes, (b) SEP measurements from SOHO at different energy channels. 


3.3. The Importance of Including All CMEs 
in Quick Succession 

[ 37 ] Our operational experience shows that it is important 
to include all CMEs occurring in quick succession (e.g.. 


within a 24 h period), especially for those CMEs which over- 
lap in their propagation paths. The inclusion of all CMEs is 
to ensure realistic representation of interactions between 
CMEs and the ambient solar wind and interactions among 
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Figure 9. An example showing that the combined WSA + ENLIL + Cone simulation results 
can be used for SEP monitoring/forecasting. 


CMEs themselves. An accurate capture of aforemen- 
tioned interactions is important for forecasting CME arrivals 
and impacts. 

[ 38 ] Within roughly a 24 h period during 18 January 2012 
and 19 January 2012, three CMEs occurred. Details are 
listed in Table 3 in a sequential order. Since CME No. 2 is 
slow (372 km/s) and narrow (half-cone angle of 27° only), 
its impact is negligible in comparison with the other two 
(at a speed of 520 km/s and 1020 km/s, respectively). 
However, it is still included in the simulation results 
shown below. 

[ 39 ] Two runs were carried out during our space weather 
operation, with Run 1 including the two CMEs (No. 2 and 
No. 3) on 19 January 2012 and Run 2 including all three 
CMEs including the one on 18 January 2012 (No. 1, No. 2, 
and No. 3). Figure 10 shows the snapshots (at 20 January 
2012T06:00Z) of the two WSA+ENLIL+Cone runs while 
still at the relatively early stage of the CMEs' evolution/ 
propagation, while Figure 11 shows the snapshots of 
the two runs when the CMEs were at the vicinity of Earth 
(at 22 January 2012T00:00Z). Figures 10 and 11 bear the 
same format as Figures 2 and 9, illustrating the rather 
standard density distribution in three different cut planes. 


[ 40 ] The predicted arrival time at Earth from Run 1 was 
21 January 2012T22:29Z and the predicated arrival time 
from Run 2 was 22 January 2012T00:52Z. By including the 
CME occurring on 18 January, Run 2 improves the fore- 
casted arrival time by 2h 23 min. Figure 12 shows the 
simulated solar wind density/speed (top panel, with the 
density in black and speed in red) and the total magnetic 
field (bottom panel) from the two runs (Figures 12a and 
12b), with the arrival time being marked by the vertical 
lines. Note that because the two runs were performed 
during our real-time space weather operations, the time 
window for two runs is slightly different. The actual 
observed shock arrival time at ACE spacecraft is 22 
January 2012T05:14Z as indicated by Figure 13. The top 
panel of Figure 13 shows the magnitude of IMF at ACE in 
black, the z component of IMF in red, and the y component 
of IMF in blue. The bottom panel of Figure 13 is the 
measured solar wind speed. The green vertical line 
indicates interplanetary shock arrival due to the CMEs. 

[ 41 ] From the simulation results of the two runs, we can see 
that the earlier and slower CME on 18 January seems to 
merge together with the faster CME on 19 January and slow 
down its arrival by 2 h and 23 min. Run 2 most likely reflects 


Table 3. The Three CMEs Occurred in Quick Succession During 18-19 January 2012 


CME Starting Time 

Speed km/s 

Direction LON/LAT in HEEQ 

Half-Cone Angle (Degrees) 

Time at 21.5 Rs 

1. 2012-01 -18T 13:25Z 

520 

3/ -24 

32 

2012-01-18T20:20Z 

2. 2012-01 -19T10:12Z 

372 

23/48 

27 

2012-01-19T20:51Z 

3. 2012-01 -19T15:10Z 

1020 

-21/46 

69 

2012-01-19T18:34 
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(a) The run with 2 CMEs 
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(b) The run with 3 CMEs 
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Figure 10. The simulated 3-D solar wind density map at 20 January 2012 06:00 UT for Run 1 (a) 
that includes two CMEs on 19 January 2012 and Run 2 (b) that includes CMEs in Run 1 and 
another CME on 18 January 2012. 
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(a) The run with 2 CMEs 
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(b) The run with 3 CMEs 
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Figure 11. Similar to Figure 10, but showing the solar wind density map at 22 January 2012 
00:00 UT for Run 1 and Run 2. 
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(a) The run with 2 CMEs 




Figure 12. The simulated solar wind density/ speed (top panel, with the density in black 
and speed in red) and the total magnetic field (bottom panel) from Run 1 and Run 2 
(Figures 12a and 12b), with the arrival time being marked by the vertical lines. The 
predicted arrival time at Earth from Run 1 was 21 January 2012 22:29 UT and that from 
Run 2 was 22 January 2012 00:52 UT. 


the actual ambient solar wind conditions more precisely. The 
interactions among CMEs with the ambient solar wind and 
other solar wind structures such as CIRs have been actively 
pursued as research topics [e.g., Gopalswamy et al, 2009]. 
Here we want to emphasize their importance in space 
weather operations, as indicated by this example (and many 
other examples not shown). As a general practice, all CMEs 
occur within a short period (roughly 24 h depending on indi- 
vidual situations) are included in the WSA+ENLIL+Cone 
simulation when assessing their arrival and impacts. 

3.4. The Importance of Other Parameters 

[ 42 ] In addition to those mentioned above, a number of 
other important factors must be taken into consideration. 

[ 43 ] An accurate prediction of CME propagation and 
transport in interplanetary space requires an accurate 
portrait of CME parameters and the medium through 
which it travels. This paper's intent is not to dismiss previ- 
ous research efforts in addressing these two essential 
aspects of the CME forecasting, but rather to reiterate their 
importance and to urge further efforts by the community to 
improve forecasting capability. 

[ 44 ] Among a CME's kinetic attributes, its propagation 
direction has particular importance. It also tends to possess 
large uncertainties (for example, in comparison to speed 
and angular width) when determined operationally. The 
importance of CMEs' propagation direction to an object 


of impact in interplanetary space can be readily under- 
stood. Research results of some time ago [e.g., De Young 
and Hundhausen, 1973; D'Uston et al., 1981] demonstrate that 
CMEs are highly directional and that their effects at 1 AU 
differ substantially depending on their propagation direc- 
tion. However, the terminology of that era (everything 
was referred to as "flare" associated, and no distinction 
was made between flares and CMEs) limits their recog- 
nition as studies of CMEs. In addition, for people who are 
involved in operations, obtaining accurate information of 
the CME propagation direction can present practical 
concerns. Uncertainties in CME propagation direction are 
almost inevitable as data gaps in real-time coronagraph 
images are common and concurrent images from more 
than one vantage point are sometimes not possible. The 
CME direction not only determines how much an object 
of interest is affected and whether it is affected at all (flank 
impacts or a miss), but also determines what kind of ambi- 
ent solar wind through which the CME propagates. So the 
influence of the propagation on the CME transport is evi- 
dent. Simulation results of six WSA+ENLIL+Cone runs in 
Table 4 serve as a good example. Roughly based on the 
CME that occurred on 11 April 2013, an ensemble of runs 
was conducted to examine how the propagation affects its 
arrival time at Earth. The input parameters of the CME 
for each run are tabulated in Table 4. For any two runs with 
the same speed (1000 km/s, or 900km/s, or 800km/s), 10° 
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Figure 13. The top panel shows the magnitude of IMF at ACE in black, the z component of 
IMF in red, and the y component of IMF in blue. The bottom panel is the measured solar wind 
speed. The interplanetary shock arrival due to the CMEs (at 22 January 2012 05:14 UT) is 
indicated by the green vertical line. 


longitudinal difference results in roughly 2 h difference in 
arrival time. It should be mentioned that the modeled 
CME for the comparative analysis is either 45° or 35° east 
off the Sun-Earth line. If the simulated CME has a more di- 
rect hit on Earth, the difference in arrival time could be big- 
ger. Our results shown here serve as further confirmation 
of previous results [e.g., De Young and Hundhausen, 1973; 
D'Uston et al., 1981]. 

[ 45 ] Table 4 also indicates that CME speed is an impor- 
tant factor in determining its arrival time. 

[ 46 ] Since CMEs have to propagate through the ambient 
solar wind, an accurate description of the ambient solar 
wind is also crucial. Even though this topic has been in- 
vestigated extensively before [e.g.. Case et al., 2008 and 


references therein], how to best represent the ambient solar 
wind in an operational setting remains a challenging task. 

[ 47 ] A detailed investigation on how various CME param- 
eters affect the arrival time will be left for a future study. 

4. Summary and Discussion 

[ 48 ] Better forecasting of a CME's propagation/ evolution 
boils down to an improved representation of properties of 
the CME itself and the medium it travels through. Based 
on our operational experience, we have identified a few 
important considerations worthy of sharing, not only with 
regard to modeling CMEs, but also concerning the use of 
CME modeling results for additional benefits. Even though 


Table 4. Effects of Propagation Direction on the Arrival Time of a CME 


Time @ 21.5 Rs 

Speed (km/s) 

Direction LON/LAT in HEEQ 

Half-Cone Angle (Degrees) 

Arrival Time @ 1 AU 

2013-04-llT10:50Z 

1000 

-35/0 

55 

2013-04-13T08:19Z 

2013-04-llT10:50Z 

1000 

-45/0 

55 

2013-04-13T10:12Z 

2013-04-llT10:58Z 

900 

-35/0 

55 

2013-04-13Tll:45Z 

2013-04-llT10:58Z 

900 

-45/0 

55 

2013-04-13T13:37Z 

2013-04-llTll:03Z 

800 

-35/0 

55 

2013-04-13T15:38Z 

2013-04-llTll:03Z 

800 

-45/0 

55 

2013-04-13T17:26Z 
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existing observations relevant to CMEs are not ideal (and 
more is still needed), they are much better than before with 
the measurements made available from the recent launch 
of STEREO A, STEREO B, and SDO. Maximizing all the 
available information to get the best assessment of CME 
properties is emphasized. Besides the conventional use of 
WSA+ENLIL+Cone model results for CME arrival/impact 
forecasts, the simulation results can also provide a visual 
representation of the distribution of SEPs, which is stressed 
in the paper. The third important point is to include all 
CMEs occurring in quick succession where the time win- 
dow for this guideline is set roughly as 24 h. In addition, 
we have also identified a number of factors that affect 
CME arrival forecasting with a particular emphasis on the 
CME propagation direction. 

[ 49 ] The coupled WSA+ENLIL+Cone model suite, as 
a sophisticated, physics-based 3-D model, is not only 
capable of performing tailored, complex calculations for 
research purpose, but also proves itself a powerful tool 
for space weather operations. Comparing to other global 
3-D models, its fast running time (about 20 min each run) 
is another advantage. As mentioned above, however, it also 
has its own limitations. For example, a CME is modeled 
as an overpressured plasma cloud without its internal 
magnetic structure(s); shocks from CMEs are modeled to 
a limited degree. 

[ 50 ] Real-time modeling and/or forecasting SEPs accu- 
rately continue to present one of the most challenging tasks 
in space weather research. While we should not ignore the 
added benefit of using the ICME modeling results to qual- 
itatively assess where SEPs are expected in interplanetary 
space, caution needs to be exercised regarding its limita- 
tions. There are very uncertain relationships between 
shocks and SEPs as observed in many cases. The particle 
intensity and energy level of SEPs can vary greatly at 
shocks of nearly identical Mach number. The shock 
geometry, angle between B and the shock normal, and 
the required seed population of superthermal particles all 
play important roles in SEP acceleration [e.g., Reames, 
1999, Tylka and Lee, 2006; Tylka et al., 2006; Li et al., 2009]. 
Shock strength alone cannot provide a very accurate pre- 
diction of SEPs. Thus, the MHD model shock strength only 
allows for a rough estimate of where SEPs might occur. If 
coupled with a particle transport code such as SEPMOD 
[Luhmann et al., 2007; Mewaldt et al., 2013], the extended 
WSA+ENLIL+Cone model could improve the SEP model- 
ing capability. However, such SEP modeling results are 
not yet routinely available. 

[ 51 ] To account for uncertainties in the initial CME 
parameters, an ensemble approach has been introduced 
and explored by our team members as a possible future 
direction [Pulkkinen et al., 2011]. It holds the promise 
of meaningfully characterizing the uncertainty involved 
in the CME modeling/forecasting process by providing 
a dynamic error bar. It should be noted, however, that 
accurate determination of the mean values of the CME pa- 
rameters remains crucial in the ensemble approach. Until 
continuous operational data streams for CME parameter 


determination become a reality, scientific expertise within 
a forecasting center/team (and constant knowledge trans- 
fer from members who have the expertise to other team 
members) always needs to be invoked in forecasting 
CMEs. 
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